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Abstract There are many localities yielding the Pleistocene mammalian faunas in China. It 
offers excellent material for the study of mammalian evolution, biochronology, paleoecology, 
paleoenvironment, paleozoogeography, etc. Faunal assemblage characters and taxonomic 
extinction rates were widely used for determining the faunal ages in biochronology. Faunal binary 
similarity coefficients sequenced according to Brainerd-Robinson’s rule and antiquity coefficients 
were further developed methods in biochronology for dating the faunal ages. The faunal binary 
similarity coefficients are based on the presence or absence of a taxon in a fauna. It is similar to 
the presence or absence of a character of a species in cladistic analyses for phylogeny, and all 
faunas have a special ancestor-descendant relationship. The present work is an attempt to find the 
relationship of the faunas with eladistie,methods by selecting three groups of faunas sequenced 
by faunal binary similarity coefficients/aceérding to Brainerd-Robinson’s rule and antiquity 
coefficients, to compare the results with different/methods, and then to estimate the ages of the 
faunas not yet dated by physical or chemical methods..The estimations are as follow: Gulongshan 
in Dalian, Liaoning Province: 16—20 ka; Shanchengzi at Benxi, Liaoning Province: 20-30 ka; 
Hualongdong at Dongzhi Man site, Anhui Province: 150—400 ka; Xinghuashan at Nanzhao Man 
site, Henan Province: 150-400 ka; Donghe at Luonan Man site, Shaanxi Province: 500—700 
ka; Bailongdong at Yunxi Man site, Hubei Province: 500-850 ka; Meipu at ‘Yunxian, Hubei 
Province: 500-850 ka; Mohui at Tiandong Man site, Guangxi Autonomous Region: 1.28M a4 
Juyuandong at Liucheng, Guangxi Autonomous Region: 1.2—1.5 Ma; Chutoulang at Chifeng, Nei 
Mongol Autonomous Region: 1.6-1.9 Ma; Renzidong at Fanchang, Anhui Province: 1.9-2.4 Ma. 
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1 Introduction 


Cladistics is a method of hierarchical analysis in which a branching diagram or 
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cladogram is used to group entities as nested sets based on special common similarities (Young, 
1986). It is widely used for searching phylogenetic relationship of given groups of organisms 
(Hennig, 1965; Ashlock, 1974; Nelson and Platnick, 1991; de Pinna, 1991). It is also used 
in other disciplines other than phylogeny of organisms, such as reconstruction of ancient 
continents evolution (Young, 1986), determination of the ages of divergent nodes of different 
clades (Ashlock, 1974). In short, any group of individuals that can be hypothesized to have 
a common ancestor, and to which a set of shared characteristics may or may not apply, can 
be compared in pairs (Platnick and Cameron, 1977). Thus, cladistics is also used in textual, 
linguistic, comparative mythology and folktale (Tehrani, 2013; d'Huy, 2013; Ross et al., 2013), 
stemmatology (Robinson and O’Hara, 1996; Kapitan, 2017), and even in astrophysics (Fraix- 
Burnet et al., 2006). 

The Pleistocene mammalian faunas are very numerous in China and they offer very 
usefulmaterial for the study of mammalian evolution, biochronology, paleoecology, 
paleoenvironmentmpaleozoogeography, and so on. Some faunas have been dated with physical 
or chemical methods*f6m, absolute ages, but some others have not yet and in need of dating. 
During the study of the’mammalian fauna from Xiaogushan site (Dong et al., 2010), the faunal 
binary similarity coefficients sequenced according to Brainerd-Robinson’s rule, a method 
proposed by Chen (1983, 2005), wasAntroduted to find the relationship between Xiaogushan 
fauna and some selected faunas, especially fombiochronological sequence of the faunas. 
The method was developed further into antiquity coefficients to find the biochronological 
sequence of the faunas during the establishment of*biochronological framework of Homo 
erectus horizons in China (Dong, 2016) and the study of the Early Pleistocene fauna from 
Chutoulang in northeastern China (Dong et al., 2017). The faunal binaty similarity coefficients 
are based on the presence or absence of a taxon in a fauna. It is similar to/the presence or 
absence of a character of a species in cladistic analyses for phylogeny in mathematiessense, 
and all faunas have their special ancestor-descendant relationship. It is therefore possible 
to use cladistic methods to find the affinity of the faunas and to interpret biochronological 
sequences, zoogeographical boundaries and paleoenvironmental evolution of the faunas. Here 
we show a test of application of cladistic methods to find the relationship of three groups of the 
Pleistocene faunas and to learn the effectiveness of the methods, and give our estimated ages 
on the faunas not yet dated by physical or chemical methods. 


2 Material and methods 


Three groups of the Pleistocene mammalian faunas analyzed by binary similarity 
coefficients in above mentioned three studies (Dong et al., 2010, 2017; Dong, 2016) were 
selected here again for cladistic analyses. 

The first group includes Xiaogushan fauna (Dong et al., 2010) and nine faunas with 
comparable geographic and chronological significances, e.g. Gulongshan fauna (Zhou et al., 


201906.00028v1 


chinaXiv 


ChinaXiv@ F RAT! 


Dong et al. — Cladistic approach on chronological relationship of the Pleistocene mammalian faunas 3 


1990) from Dalian, Liaoning; Miaohoushan fauna (Zhang et al., 1986) and Shanchengzi fauna 
(Zhang et al., 1986) from Benxi, Liaoning; Yushu fauna (IVPP, 1959) from Yushu area, Jilin; 
Yanjiagang fauna (NET et al., 1987) from Haerbin, Heilongjiang; Dingcun fauna (Pei, 1958) 
from Xiangfen, Shanxi; Upper Cave fauna (Pei, 1940) from Zhoukoudian, Beijing; Lingjing 
fauna (Li and Dong, 2007; Zhang et al., 2009) from Xuchang, Henan; and Salawusu (Sjara- 
osso-gol) fauna (Qi, 1975; Tong et al., 2008) from Erdos, Nei Mongol. The ages of the faunas 
in this group range in the Late Pleistocene except that of Miaohoushan which is of the late 
stage of the Middle Pleistocene. The selected faunas total 10 and referred taxa 116 (Appendix I). 

The second group are faunas from Homo erectus horizons in China, e.g. Danawu 
horizon from Yuanmou Man site (Lin et al., 1978), Mohui horizon from Tiandong (Wang et 
al., 2007), Gongwangling horizon from Lantian Man site (Hu and Qi, 1978), Longgudong 
horizon from Jianshi Man site (Zheng, 2004), Quyuanhekou horizon from Yunxian Man site 
(Ji, 999)», Meipu horizon from Yunxian (Xu, 1978), Donghe horizon from Luonan (Xue 
et al., 1999), Bailongdong horizon from Yunxi Man site (Wu et al., 2009), Longtandong 
horizon from HexianyMan site (Huang et al., 2012), Chenjiawo horizon from Lantian (Qi, 
1989), Qizianshan horizon from Yiyuan Man site (Lü et al., 1989), Xinghuashan horizon from 
Nanzhao Man site (Qiu et al.,1982), Huludong horizon from Nanjing Man site (Wu et al., 
2002) and Hualongdong horizon from Dongzhi Man site (Gong et al., 2014). Since there are 9 
layers yielding Homo erectus material (Li andi, 1981) and the hominin horizons at Locality 
1 of Zhoukoudian Site were defined according tosthree cycles of sedimentations related to 
environmental changes (Xu et al., 1999), Layers*8—-11, dzayers 6-7 and Layers 3—5 were 
considered as lower, middle and upper Homo erectus horizons tespectively. The ages of the 
faunas in this group range from the middle stage of the Early Pleistocéne/to the middle stage 
of the Middle Pleistocene. Because some faunas were collected with screen washing technique 
for small mammals, but the others not, and in order to reach balanced results wejusteselect the 
taxa of large mammals from these faunas. The faunas from all Homo erectus horizons from 
China total 17 and referred taxa 123 (Appendix II). 

The third group are all Early Pleistocene faunas with the number of the taxa above 
30 from northeastern, northern, northwestern, central and southern parts of China such as 
Chutoulang (Zhang, 1989; You and Zhang, 1989; Dong et al., 2017), Nihewan (s.s.) (Teilhard 
de Chardin and Piveteau, 1930; Qiu, 2000; Tong et al., 2011), Gongwangling (Hu and Qi, 
1978), Longdan (Qiu et al., 2004), Renzidong (Jin and Liu, 2009), Longgupo (Huang et al., 
1991), Longgudong (Zheng, 2004), Juyuandong (Gigantopithecus Cave) (Pei, 1987; Han, 
1987; Han and Xu, 1989) and Dadong (Jin et al., 2009). A French fauna from upper Saint- 
Vallier (LD3) about 1.5 km northeast of Saint-Vallier Village (45°11'16.94"N, 4°49'47.69"E, 
360 m above sea level) (Guérin et al., 2004) was also selected as a classic European 
Villafranchian fauna for comparison. The ages of the faunas in this group range from the early 
stage of the Early Pleistocene to the middle stage of the Early Pleistocene. The selected faunas 
total 10 and referred taxa 317, small mammals included (Appendix III). 
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All three groups were then gathered together as a single assembled group ranging from 
the early stage of the Early Pleistocene to the late stage of the Late Pleistocene for cladistic 
analyses. Since the faunas from Gongwangling and Longgudong are present in both the second 
and third groups, the faunas in assembled group total 35 and referred taxa 418, small mammals 
included (Appendix IV). 

All extinct taxa, including regional extinct ones, were considered as a presumed archaic 
fauna and those of extant taxa as a presumed modern fauna in every group. 

The presence or absence of a species, including conformis species, in a fauna is 
considered as equivalent to the presence or absence of a character in a taxon in phylogenetic 
cladistic analyses and marked by “1” (present) or “0” (absent) in the data matrices in the 
present work. The question mark “?” in the data matrix indicates that the taxon of the fauna 
is/likely.close to the determined species, and it is equivalent to the missing character in 
phylogenetic eladistic analyses. 

The data niatrices (Appendices I-IV) were edited with MS Excel 2016 and NEXUS Data 
Editor (Version 0.4.8 by Roderic D. M. Page) and then processed with PAUP* (Version 4.0a 
build 164 for 32-bit Microsoft Windows by David L. Swofford). The presumed archaic fauna 
was set as outgroup for the analyses, and the presumed modern fauna was considered as a 
crown group. All “characters” are of unordered type with equal weight. Optimality criterion 


was set as “parsimony” with heuristic search./The optimal trees and consensus trees were 


Md Sayed and then organized with TreeViewX 
GLS (Version 05.Qyby Roderic D. M. Page). 
SDD Thé results of cladistic analyses were 


SCZ then compared with those ofbinary faunal 
XGS similarity coefficients sequeneed, according 


SLWS to Brainerd-Robinson’s rule (Chen, 1983, 

DC 2005) and antiquity coefficients (Dongs 

YJG 2016), extinction rates, as well as dated ages, 

YS to test the effectiveness of the analyses on 

Si faunas and estimate the ages of the faunas 

ii not yet dated by any physical or chemical 
rc 


Fig. 1 Optimal tree for selected Late Pleistocene 
mammalian faunas from China 
Arc. Presumed archaic fauna; DC. Dingcun at 
Xiangfen; GLS. Gulongshan in Dalian; LJ. Lingjing in 
Xuchang; Md. Presumed modern fauna; 

MHS. Miaohoushan at Benxi; SCZ. Shanchengzi at 
Benxi; SDD. Upper Cave at Zhoukoudian; 
SLWS. Salawusu (Sjara-osso-gol) in Ordos; 
XGS. Xiaogushan at Haicheng; YJG. Yanjiagang in 
Harbin; YS. Yushu at Yushu area 


methods. 


3 Results of cladistic analyses 


For the first group of the faunas, the 
heuristic tree search tried 902 rearrange- 
ments and found 244 best trees in which one 
optimal tree with 10 nodes was retained, in 
which 9 nodes have dichotomous branches 
and only one has three (Fig. 1). 
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For the second group of the faunas, 
the heuristic tree search tried 48498 
rearrangements and found 275 best trees in 
which 10 optimal trees were retained. The 
strict and semi-strict consensus trees of the 
retained optimal trees are the same. Both 
of them have 11 nodes with dichotomous 
branches, one node with three branches 
and another with four (Fig. 2, right). The 
Adams consensus tree has some minor 
differences from the strict and semi-strict 
consensus trees. It has two nodes with 
thre@*branches and the rest 14 nodes have 
dichotomous branches (Fig. 2, left). 

For the third group of the faunas, the 
heuristic tree search“tried 3737 rearrange- 
ments and found 620 best treg§ ña which 
seven optimal trees were retained. The 


CTL CTL 
NHW NHW 
LD Md 
StV StV 
Gwl LD 
DD DD 
JYD Gwl 
Md JYD 
Lgd Lgd 
LGP LGP 
RZD RZD 
Arc Arc 


Fig. 3 Adams (left) and strict consensus (right) 
trees for selected Early Pleistocene mammalian 
faunas from China 
Arc. Presumed archaic fauna; CTL. Chutoulang 
in Chifeng; DD. Dadong at Chongzuo; 

Gwl. Gongwangling at Lantian; JYD. Juyuandong 
(Gigantopithecus Cave) at Liucheng; 

LD. Longdan in Linxia Basin; Lgd. Longgudong 
at Jianshi; LGP. Longgupo at Wushan; 

Md. Presumed modern fauna; NHW. Xiashagou 
in Nihewan Basin; RZD. Renzidong at Fanchang; 
StV. Saint Vallier in France for comparison 


Md Md 

Ltd Ltd 
Cjw Cjw 
Zkd 3-5 Zkd 3-5 
Zkd 8-11 Zkd 8-11 
Zkd 6-7 Zkd 6-7 
Qzas Qzas 
Ts Ts 

Xhs Xhs 

Dh Dh 
Qyhk Hid 
Gwl YmM 
Hid Qyhk 
YmM Gwl 

Mp Mp 

Bid Bld 
Lgd Lgd 

Mh Mh 

Arc Arc 


Fig. 2 Adams (left) and strict consensus (right) trees for 
the faunas from Homo erectus horizons in China 
Arc. Presumed archaic fauna; Bld. Yunxi Man locality at 
Bailongdong; Cjw. Lantian Man locality at Chenjiawo; 
Dh. Luonan Man locality at Donghe; Gwl. Lantian Man 
locality at Gongwangling; Hld. Dongzhi Man locality at 
Hualongdong; Lgd. Jianshi Man locality at Longgudong; 
Ltd. Hexian Man locality at Longtandong; Md. Presumed 
modernefauna;¢Mh. Tiandong Man locality at Mohui 
Cave; Mp. Yunxian Man locality at Meipu; 

Qyhk. Yunxian Man locality at Quyuanhekou; Qzas. 
Yiyuan Man localities at Qizianshan; Ts. Nanjing Man 
locality at Huludong; Xhs. Nanzhao Man locality at 
Xinghuashan; YmM. Yuanmou Man localitysat Danawu; 
Zkd 3-5, Zkd 6-7, Zkd 8-11. Layers 3—5{ 637, 8-11 at 
Locality 1 of Zhoukoudian Peking Man Site 


strict and semi-strict consensus trees of the seven 
optimal trees retained are the same, each tree has 
only three nodes: the first one has three branches, 
the second one has 9 branches and only the last 
one has dichotomous branches (Fig. 3, right). 
The Adams consensus tree is more distinct. 
It has seven nodes, in which four nodes have 
dichotomous branches, two nodes have three 
branches and one has four (Fig. 3, left). 

For the assembled group of the faunas, the 
heuristic tree search tried 1452544 rearrange- 
ments and found 995 best trees in which 169 


optimal trees were retained. The strict and semi- 
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strict consensus trees of 169 trees are very similar to each other; the former has one more 
node with five branches compared with the latter. The semi-strict consensus tree has 25 nodes, 
in which 20 nodes have dichotomous branches, one has three, two have four and the other 
two have five (Fig. 4, right). The Adams consensus tree has 29 nodes, in which 24 nodes 
have dichotomous branches, three nodes have three branches and two have four (Fig. 4, left). 
Miaohoushan (MHS) fauna of the first group appears in the clade formed by many faunas from 
the second group, and some faunas from the second group, e.g. Mohui (Mh) and Longgudong 
(Lgd), appear on the stem formed by many faunas from the third group. 


Md Md 
GLS GLS 
SDD SDD 
SCZ SCZ 
XGS XGS 
SLWS YJG 
YJG YS 
YS SLWS 
pc DC 

LJ LJ 
Zkd 3-5 Zkd 3-5 
Zkd 8-11 Zkd 8-11 
Zkd 6-7 Zkd 6-7 
Ltd Ltd 
MHS MHS 
Qzas Qzas 
Ts Ts 
Xhs Qyhk 
Qyhk Gwl 
Gwl Xhs 
Cjw Cjw 
Dh Dh 
Hid CTL 
CTL NHW 
NHW LD 
LD StV 
StV hiid 
Mp Mp 
YmM YmM 
Bid Bld 
JYD JYD 
DD DD 
Mh Mh 
LGP LGP 
Lgd Lgd 
RZD RZD 
Arc Arc 


Fig. 4 Adams (left) and semi-strict consensus (right) trees 
for selected Pleistocene mammalian faunas from China 
See Figs. 1-3 for abbreviations of the localities 


4 Discussion 


4.1 The first group of faunas 

The result of cladistic analysis on the first group of the Pleistocene faunas is 
chronologically and geographically significant (Fig. 1; Table 1). The branches of the crown 
clade such as Xiaogushan (XGS), Shanchengzi (SCZ), Upper cave (SDD), Gulongshan (GLS) 
are geographically close to each other and chronologically related (Fig. 1). It is the same for 


201906.00028v1 


chinaXiv 


ChinaXiv@ F RAT! 


Dong et al. — Cladistic approach on chronological relationship of the Pleistocene mammalian faunas 7 


the clade by Yanjiagang (YJG) and Yushu (YS). For the clade by Dingcun (DC) and Salawusu 
(SLWS), they are separated by a moderate distance of about 360 km, but chronologically 
very close to each other (Table 1). The result is very close to that by binary faunal similarity 
coefficients sequenced according to Brainerd-Robinson’s rule (Dong et al., 2010). It is also 
close to that by antiquity coefficients and extinction rates (Table 1). The cladistic analysis 
appears effective for this group of the faunas. Gulongshan and Shanchengzi faunas have not 
yet dated by any physical or chemical methods. Furthermore, the result of the analyses of the 
first group is also very close to that with assembled group. Based on the cladistic analyses, 
Gulongshan fauna is closest to presumed modern faunas (Figs. 1, 4), it is the same by binary 
faunal similarity coefficients sequenced according to Brainerd-Robinson’s rule (Dong et al., 
2010). But it is slightly older than Upper Cave fauna by antiquity coefficients and extinction 
tates. (Table 1). Shanchengzi fauna is placed on the stem between the branches of Xiaogushan 
and(UppemCave faunas by cladistic analyses (Fig. 1), in accordance with the results by binary 
faunal similarity,coefficients sequenced according to Brainerd-Robinson’s rule and antiquity 
coefficients (Deng @tal., 2010). Miaohoushan (MHS) fauna is close to the presumed archaic 
fauna of the first group (Fig. 1) or close to the Middle Pleistocene faunas (Fig. 4), and it is 
the oldest fauna in the first grotip\with all methods (Table 1). Lingjing (LJ) fauna appears as 
the second oldest fauna in the first/group/with all methods (Fig. 1; Table 1). For the faunas 
without numeric dating, we estimate the age of Gulongshan (GLS) fauna roughly equivalent 
to that of Upper Cave one, as 16—20 ka. We estimated,the age of Shanchengzi fauna as 20-30 
ka, evidently older than that of Upper Cave fauna, overlaps that of Yanjiagang fauna, partly 
overlaps that of Yushu fauna and younger than that of Xiaogushan fauna (Table 1), although 
the extinction rates indicate its age would be between those of Yushu and Xiaogushan. 


Table 1 Comparison of the results of the 1“ group faunas with different methods 


Cladistics B.-R.’s rule? Antiquity” Extinction rate” Dated age Estimated,age 
GLS GLS SDD (100) SDD (13.64%) GDS; 16-20’ka 
SDD SDD GLS (120) GLS (17.65%) SDD: 18 ka (Huang, 1991) 

SCZ SCZ SCZ (126)  YJG (26.67%) SCZ: 20-30 ka 
XGS XGS SLWS (128) SLWS (27.27%) YJG: (27.0+0.6) ka (Wei et al., 1986) 
SLWS YS YJG (132) DC (30.00%) YS: 27-40 ka (Sun et al., 1981) 

DC YJG DC (144) YS (33.33%) XGS: 30-80 ka (Zhang et al., 2010) 

YJG SLWS YS (146) SCZ (35.9%) | SLWS: 61—68 ka (Yin and Huang, 2004) 

YS DC XGS (156)  XGS (37.50%) DC: 100 ka (Huang, 1991) 

LJ LJ LJ (162) LJ (52.94%) LJ: 105-125 ka (Li et al., 2017) 

MHS MHS MHS (174) MHS (53.13%) MHS: <400 ka (NET et al., 1987) 


B.-R. Brainerd-Robinson (Chen, 1983). 1) Dong et al., 2010. See Fig. 1 for abbreviations for fauna localities. 


4.2 The second group of faunas 


The results of the cladistic analyses for the second group (Fig. 2) show that Hexian Man 
fauna (Ltd) is closest to the presumed modern fauna of the group and they form a clade with 
Chenjiawo fauna. Three faunas of Peking Man horizons at Zhoukoudian form a reasonable 
clade, but surprisingly the faunas from upper and lower horizons are more closely related than 
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the middle one, the result is the same as that according to B.-R.’s rule (Dong, 2016). They 
form a clade successively with the faunas from Yiyuan Man site (Qzas), and then Nanjing 
Man site (Ts). The faunas from Gongwangling (Gwl) and Yunxian Man site (Qyhk) also form 
a reasonable clade. But the clade by the faunas from Yuanmou Man site (YmM) and that of 
Dongzhi Man site (HId) is unexpected. It results from the numerous undeterminable species 
in Hld fauna whose genera are the same with those of YmM, and the taxon number of Hld 
is not sufficient on the other hand. Therefore, it does not indicate their close relationship 
in biochronology. It is relatively complicated in comparing the cladistic results with those 
by binary faunal similarity coefficients sequenced according to Brainerd-Robinson’s rule, 
antiquity coefficients and extinction rates (Table 2). It likely results from the zoogeographic 
differences of the faunas that make the fauna compositions very varied. We gave our estimated 
agesof the faunas not yet dated by physical or chemical methods in Table 2 based on the 
results by different biochronological methods. 


Table 2/%Comparison of the results of the 2™ group faunas with different methods 


Cladistics B.-R.’s rule” Antiquity,” Extinction rate” Dated age Estimated age 
Ltd Gwl Hid (81) Mp (33.33%) Hld: 150-400 ka 
Zkd 3-5: 230-310 ka (Zhao et al., 
1985);>400 ka (Shen et al., 2001) 
Zkd 3-5 Mh Xhs (88) Xhs (40,00%) Xhs: 150-400 ka 
Ts: 350 ka (Chen et al., 1996); 


Ciw Ltd Ts (87) Hd (36.36%) 


Zkd 8-11 Lgd Zkd 3-5 (89) Ts (43.75%) 
500 ka (Wang et al., 2002) 
Ltd: 150-190 ka (Chen et al., 1987); 
Zkd 6-7 Bld Qzas (90) Bld (46.43%) 
300ka (Huang, P H et al., 1995) 
Qzas Qyhk Cjw (92) Zkd 3-5 (47.06%) Qzas: 320#420/ka (Han et al., 2016) 
Ts Mp Zkd 6-7 (94) Cjw (50.00%) Zkd 6-7: 380 ka (Zhao etal. 1985) 
Xhs Dh Ltd (95) Qzas (53.85%) Cjw: 650 ka (An and Ho, 1989) 
Dh Xhs Bld (97) Dh (57.89%) Dh: 500-700 ka 
Zkd 8-11: 462 ka (Zhao et al., 1985); 
Qyhk Zkd 8-11 Dh (98) Ltd (58.06%) 
750 ka (Shen et al., 2009) 
Gwl Zkd 3-5 Mp (99) Zkd 6-7 (59.09%) Bld: 500-850 ka 
Mp Zkd 6-7 Zkd 8-11 (102) Zkd 8-11 (55.26%) Mp: 500-850’ka 
Hd Ts Qyhk (104) Qyhk (65.22%) Qyhk: 936 ka (de Lumley et al., 2008) 
Gwl: 1.15 Ma (An and Ho, 1989); 
YmM Qzas Gwl (106) Lgd (72.50%) 
1.7 Ma (Zhu et al., 2015) 
Bld Cjw Lgd (114) Gwl (75.00%) Lgd: 2.14 Ma (Gao and Cheng, 2004) 
Lgd Hld Mh (117) Mh (80.00%) Mh: 1.2-1.8 Ma 
YmM: 700 ka (Hyodo et al., 2002); 
Mh YmM YmM (120) YmM (80.77%) 


1.7 Ma (Zhu et al., 2008) 


B.-R. Brainerd-Robinson (Chen, 1983). 1) Dong, 2016. See Fig. 2 for abbreviations for fauna localities or horizons. 


4.3 The third group of faunas 


The strict and semi-strict consensus trees illustrate that Nihewan (NHW) and Chutoulang 
(CTL) faunas are the closest, Renzidong (RZD) fauna is the basal fauna in the third group, and 
the other faunas are generally very close to each other (Fig. 3, right). Dadong (DD) is closely 
related to Gigantopithecus Cave (JYD), Longgupo is closely related to Longgudong (Fig. 3, 
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left; Fig. 4), just as they are closely related in zoogeography. The close relationship between 
JYD and DD (dated as 1.2 Ma, see Jin et al., 2009) is also supported by antiquity sequence 
and B.R.’s rule and we estimate the age of JYD as 1.2-1.5 Ma (Table 3). The European 
Villafranchian fauna from Saint-Vallier (StV) is the closest to that of Longdan (LD). It is also 
in accordance with the sequences by antiquity coefficients and extinction rates (Table 3). It 
might indicate that Longdan fauna and Saint-Vallier fauna were in the same zoogeographic 
region (Palaeoarctic) as well as similar biochronological period. But in terms of sequenced 
binary similarity coefficients according to B.-R.’s rule (Table 3), Longdan fauna is also very 
similar to Chutoulang and Juyuandong faunas. The close relationship between Chutoulang and 
Longdan is also found in binary faunal similarity coefficients and they are both in Palaeoarctic 
region. But with JYD of Oriental region, it can be interpreted only by their biochronological 
context (Table 3). Nevertheless, NHW and CTL form a clade with LD first and then StV in the 
assembled.group (Fig. 4). The age of CTL is very likely older than DD, younger than LD and 
StV, equivalentsor slightly younger than NHW. We estimate the age of CTL as 1.6-1.9 Ma. 
RZD is the closest branch to the presumed archaic one in both the third and assembled groups 
(Figs. 3—4). Its oldes®status in the third group is also supported by B.-R.’s rule and Antiquity 
sequence but its extinction rate”s lower than those of CTL, NHW, StV and LD, we estimate its 
age as 1.9-2.4 Ma (Table 3). 


Table 3 Comparison of the results of th@3™ group faunas with different methods 


Cladistics B.-R.’s rule? Antiquity? Extinction rate” Dated age Estimated age 

NHW Gwl Gwl(198) DD (61.1%) JYD: 1.2-1.5 Ma 
CTL Lgd DD (202) Gwl (63.4%) DD: 1.2Ma (Jin et al., 2009) 

LD DD JYD (208)  Lgd (65.9%) Gwl: 1.15 Ma (An and Hos 1989); 

1.7 Ma (Zhu et al., 2015) 

StV JYD CTL (209) LGP (69.6%) CTL: 1.6-1.9 Ma 
Gwl CTL Lgd (232) JYD (71.4%) NHW: 1.7-2.2 Ma (Liu et al., 2012) 

JYD LD LD (234) RZD (72.5%) LGP: 1.12-1.96 Ma (Huang W B et al., 1995) 

DD NHW StV (238) CTL (76.7%) Lgd: 2.14 Ma (Gao and Cheng, 2004) 

Lgd StV NHW (242) NHW (85.4%) LD: 2.16-2.55 Ma (Qiu et al., 2004) 

LGP LGP LGP (258) StV (85.7%) StV: 1.95-2.58 Ma (Sen, 2004) 

RZD RZD RZD (266) LD (96.8%) RZD: 1.9-2.4 Ma 


B.-R. Brainerd-Robinson (Chen, 1983). 1) Dong et al., 2017. See Fig. 3 for abbreviations for fauna localities. 


5 Conclusion 


Cladistic analyses are feasible methods for searching relationship of faunas. It is effective 
for comparing the faunas with large number of taxa and within the same zoogeographic 
province as the first group of the faunas for biochronological sequence. With the faunas having 
insufficient number of taxa or from different zoogeographic provinces, it is recommended to 
refer the results by other biochronological methods such as antiquity coefficients and extinction 
rates for comprehensive biochronological interpretation. 

We estimate the ages of the faunas not yet dated by physical or chemical methods based 
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on the results of present cladistic analyses and refer with those of sequences according to 
Brainerd-Robinson’s rule, antiquity coefficients and extinction rates as follow: 

Gulongshan in Dalian, Liaoning Province: 16-20 ka 

Shanchengzi at Benxi, Liaoning Province: 20-30 ka 

Hualongdong at Dongzhi Man site, Anhui Province: 150-400 ka 

Xinghuashan at Nanzhao Man site, Henan Province: 150-400 ka 

Donghe at Luonan Man site, Shaanxi Province: 500-700 ka 

Bailongdong at Yunxi Man site, Hubei Province: 500-850 ka 

Meipu at Yunxian, Hubei Province: 500-850 ka 

Mohui at Tiandong Man site, Guangxi Autonomous Region: 1.2—1.8 Ma 

Juyuandong at Liucheng, Guangxi Autonomous Region: 1.2-1.5 Ma 

GChutoulang at Chifeng, Nei Mongol Autonomous Region: 1.6-1.9 Ma 

Renzidong at Fanchang, Anhui Province: 1.9-2.4 Ma. 


Acknowledgements The present work was supported by the Strategic Priority Research 
Program of Chinese Acadeiiy of Sciences, Grant No. XDB 26000000. The authors would like 
to thank reviewers for comments and suggestions to improve the manuscript. 


Appendix can be found on the website of Vertebrata PalAsiatica (http://english.ivpp.cas.cn/sp/PalAsiatica/VP _ 
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PA BB o> E ATE SL tA S Ra AT 
BEATERS HET 


a APO AXE 2 
a PERERA ESR ARKEEN, PRR GET AE Sh ASS ORS IT AGE 100044) 
2 HEBEL Wi She Barb Ibe 100044) 
G PERFEK J6 100049) 


HE. CO PRUEA ARZTL SLO, CNEL on ik. EERE, 
TESZ., THE. HE OE TRB. TEE ERE H, 
OH A MINE KR Ze FE DT EER TE. EEN EAER A 
Brainerd-Robinsonj ll Xf 3 RERI ICALE BOE TARY, ak CHE BE BH 
BRER T E FCA RE AE I TT E OEE AA 8 ER TS 
ESIE ESR, KAFT REP EARE YE EA PEE — A 
PHASE FEES PAI. KMRL EFR RAIT ERWE ERR, 
AFE T Ei AER RET RETF RRE, KEARE 
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MAE AEA TT IS AS RIETER, 2] ET ET AE ER, PASE eG 
AR HEFT EMT AE DS I AE BS) ETT ER EDT, Fe Rs Ee 
Sh We HY 16~20 ka; WL TACK NY LIM HE 20~30 ka; BBR BAGEL A De Tle E 
74150~400 ka; PPA A A AGERE 9 46 LL TE 1 50~400 ka; PEPYS RA RAER Tl oly 
#49500~700 ka; HICH PS AHER A De oh Ye 500~850 ka; ERS EL ie oh Be 
500~850 ka; J Pa FAR Gath A eT aE 1.2~1.8 Ma; J Pa ASS Acie oh Bey 1.2~1.5 
Ma; ASE rh ap week WI SNE A 1 .6~1.9 Ma; KAS AF lel SE Wy 1.9~2.4 Ma. 

AE: HE, Heit, MPLA, FRR, EWER 
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